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Abstract 

We analyze the percolation properties of certain clusters defined on configurations of the 2- 
dimensional Heisenberg model. We find that, given any direction n, in 0(3) space, the spins almost 
perpendicular to n form a percolating cluster. This result gives indications of how the model can 
avoid a previously conjectured Kosterlitz-Thouless phase transition at finite temperature T. 
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The classical Heisenberg model in 2 dimensions (2D) 
describes the behaviour of a system of classical spins with 
short range ferromagnetic interactions [jlj . The spins are 
placed at the sites of a 2D lattice. The physics of the 
model, which has been studied both through analytical 
calculations and Monte Carlo simulations, is defined by 
equations that display a continuous 0(3) symmetry and 
it is subject to the Mermin and Wagner theorem H, i.e.: 
there are no equilibrium states with broken symmetry. 

Perturbation theory (PT) indicates that the 2D 
Heisenberg model has a critical point at zero tempera- 
ture S. Moreover, from the field-theoretical point of 
view, the spin field carries a particle of non-zero mass 
to. This mass has been calculated by applying a Bethe 
ansatz technique Q and by making a partial use of PT. 
However, PT is constructed by studying small oscilla- 
tions around the trivial configuration (all spins parallel) , 
a state which obviously violates the Mermin and Wagner 
theorem. Therefore a problem raises concerning the va- 
lidity of the above-mentioned analytical results and any 
other which relies on PT. 

The model has also been studied by numerical simu- 
lations. Among other quantities (see for example 0), 
the mass m and the magnetic susceptibility \ have been 
recently measured with good precision. The mass is eval- 
uated from the exponential decay of the 2-point corre- 
lation function G(x — x'). The magnetic susceptibility 
is extracted from the sum \ = G{x — x'). Some of 
the most recent numerical calculations are the following. 
In j(],|7j m has been determined by extrapolating the re- 
sult from small lattice sizes and large temperatures to 
smaller temperatures by using finite-size scaling. In ||] 
it has been extracted from the 2-point function by us- 
ing improved actions and very high statistics. In both 
cases agreement with the analytical calculation of to ||] 
is found within 2-4%. 

Another scenario for the 2D Heisenberg model has 
been put forward in f9|-|llj|. The model would undergo 
a Kosterlitz-Thouless (KT) phase transition at a finite 
temperature Tkt and no massive particle would be car- 
ried by the spin field in the low temperature phase. This 
scenario is a mimic of what is known to happen in the 



2D 0(2) or XY model JT2J . If T KT is low enough, a nu- 
merical simulation is not able to detect it (in HpH] it was 
argued that Tkt is indeed much smaller than the typical 
temperatures utilized for thermalization in present-day 
simulations). In |14| a finite-size analysis of the helic- 
ity modulus is used to rule out such a KT transition 
for T > 0.1. In it was also shown that the data for 

the correlation length and the magnetic susceptibility for 
temperatures T > 0.53 do not scale as the KT scenario 
predicts. 

In the present paper, we want to tackle directly the 
arguments of J^,[nJ where the percolation properties of 
certain clusters are analyzed and, after assuming a set of 
hypotheses, it is concluded that the magnetic suscepti- 
bility diverges which is sufficient to prove that the mass 
to vanishes. 

We realize the classical spins of the Heisenberg model 
by 3-component scalar fields of modulus 1 placed at each 

site x of a square lattice, $ x with (<f>x^ = 1- These fields 

interact through a nearest neighbour (n.n.) coupling and 
the hamiltonian can be written as ((xy) stands for two 
generic n.n. sites) 



(i) 



The partition function at a temperature T is Z = 
} ex P {H/T). The hamiltonian (Q) and the parti- 
tion function are invariant under 0(3) rotations. 

In the following we recall the arguments of Jl^,|ll| . Let 
us consider a configuration for this model thermalized at 
a given temperature T. Let n be an arbitrary unit vec- 
tor in the internal space of the 0(3) symmetry. To any 
such a vector we can associate various types of clusters 
on the configuration. If A is one such cluster then its 
size, defined as the number of sites contained in it, shall 
be denoted by Ca- On the other hand, its perimeter, 
defined as the number of sites along the border of the 
cluster, will be called Bj^. If a set of clusters completely 
cover the whole lattice volume with no overlap then we 
say that we have a "cluster system" . 

The Fortuin-Kasteleyn clusters (hereafter called 
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T) Eq] are made of sites connected by the bonds which 
survive the deletion process performed with the proba- 
bility 

P xy = cxp ^min{0 , ($ x ■ nj ($ y ■ nj . (2) 

The average size of the T clusters satisfies (Cp) = n\ 
where x is the magnetic susceptibility and n is a constant 
k < 1. The brackets (•) indicate average over configu- 
rations or equivalently the average calculated with the 
Boltzmann weight of the partition function. In physical 
terms, T clusters are regions of correlated spins. The set 
of J- clusters form a cluster system. 

Other clusters associated to an arbitrary unit vector n 
are the Tt + , Or and S clusters. All n.n. spins on a ther- 
malized configuration at a temperature T tipically satisfy 

110a; — <fiy 1 1 < e with a parameter e of order O 

Then for any site x the scalar product (d> x ■ nj can be 

either a) (^<fi x ■ nj > e/2, and x belongs to an Ti + cluster, 

or b) (d> x ■ nj < —e/2, so x belongs to an Hr cluster, or 

c) \4> x -n\ < e/2 and consequently x belongs to an S clus- 
ter. In simple words, the S clusters are constituted by 
sites whose spins almost lie in the plane perpendicular to 
n; on the other hand the Ti + ("H~) clusters contain sites 
whose spins are almost parallel (antiparallel) to n. The 
set of H ± and S clusters form a cluster system. 

By using a variant of the 0(3) model that includes 
a Lipschitz continuity condition, \\<f> x — (fiy\\ < S for all 
n.n. x, y (this condition does not change the physical 
properties of the model as long as S > e), it is possible 
to prove that the 7i clusters lie entirely inside the T 
clusters fl^Jll|l . Moreover no 7i + cluster has a common 
frontier with any Ti~ cluster, they are always separated 
by spins belonging to S type clusters. 

Following Ref. JT^JTTJ] , three hypotheses are now nec- 
essary to prove that the magnetic susceptibility diverges. 
The first one states the impossibility to have simultane- 
ous percolation of two or more disjoint clusters in a clus- 
ter system defined on the 2D Heisenberg model. This 
fact is known to be true in several models |L6]]. The sec- 
ond hypothesis extends the validity of the Mermin and 
Wagner theorem to hamiltonians which, like the one in- 
cluding the Lipschitz condition, are not differentiable in 
the fields. These two conditions altogether prevent the 
TC clusters from percolating. 

The third hypothesis is to assume that the S type clus- 
ters do not percolate either. Consequently none of the 
H ± or S clusters percolate. When a cluster system does 
not contain a percolating cluster then at least two kinds 
of clusters have a divergent average size [|l7) . In the clus- 
ter system of 7i ± and S, this result means that either 
(C H +) = oo or (C H -) = oo or both (actually both, to 
avoid the breaking of the second hypothesis). Moreover, 
as the H ± clusters are entirely included inside the T 



clusters, these clusters must also have a divergent aver- 
age size. Therefore the magnetic susceptibility \ diverges 
and the theory is massless in clear contrast to the calcu- 
lation of %. 

The above conclusion can be avoided if at least one of 
the three hypotheses fails. In particular we have checked 
the third one finding that there does exist a S cluster 
which percolates. In the present paper we report the re- 
sults of a numerical simulation of the Heisenberg model 
at T = 0.5 on very large lattices to give evidence that 
this is indeed the case at that temperature. We describe 
how this percolation takes place in the system, showing 
that a percolating S cluster is found for every direction n 
in the symmetry group space 0(3). Moreover all (perco- 
lating or not) S and 7i ± clusters present a high degree of 
roughness. More specifically, every spin in a cluster has 
on average one n.n. lying on the border of the cluster, (on 
a 2D square lattice each spin has four n.n.). In the fol- 
lowing we describe the technical details of our simulation 
and give the results obtained. 

We have simulated the system described by the hamil- 
tonian in Eq.(^) on several lattice sizes L 2 with L = 1024, 
1250, 1550 and 2048 at a temperature T = 0.5. The mea- 
surements and cluster analysis have been performed on 
a sample of between 1000 and 2000 configurations which 
were separated by 200 overheatbath updating steps [Gl| , 
(this updating method applies a usual heatbath step to 
the dynamical variables the angle formed by the spin 
field <fi x and the sum of their n.n.- and fixes the other 
variables by maximizing the distance between the old 
and the new updated fields <fi x in order to hasten the 
decorrelation). 

The percolation of clusters is a well-defined concept on 
infinite lattices. Obvious computer limitations force us 
to work on finite volume systems. We have overcome this 
difficulty by working with a large ratio L/£ where £ is a 
typical correlation length of the system. In this way we 
expect that the would-be percolating cluster will show 
up as a very large cluster, much larger than the others. 
We have chosen our simulation parameters in such a way 
that Lji > 4.5, see Ref. 0. 

The first task is to decide the value for e to be used in 
the construction of clusters, e should not be too large in 
order to avoid a certain (but not instructive for our prob- 
lem) percolation of the S clusters, the estimate e w \/2T 
being an appropriate value. In Fig. [l] we show the prob- 
ability distribution of values of \\<fi x — <fi v \\ for all n.n. x, 
y on a lattice of size L = 1024 at T — 0.5. The error 
bars in this and in the subsequent figures are very small: 
we show only the error bar of the data point on the top 
of Fig. 0. Motivated by this distribution we ha ve taken 
c = 1.05 which satisfies the estimate e s» \/2T = 1. We 
have checked that, by using this value of e, the Ti + and 
Or clusters touch each other only in a 0.5% of all bonds 
and this is shown in the inset of Fig. [l] where we give the 
fraction of bonds which join Ti + and HT clusters, Fh, as 
a function of e. At e = 1.05 we will already find evidence 
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for the existence of percolating S clusters. For smaller 
epsilons one cannot longer prove that the Tt^ clusters lie 
entirely inside the T clusters, which is an essential ingre- 
dient for the arguments given in Jl(],|ll| . We could choose 
smaller values of e only at lower temperatures. 

We have taken several choices for the vector ft. In all 
cases we observed exactly the same results. In this paper 
we report the figures obtained with n = (1, 0, 0). 



2.0 




0.0 0.5 1.0 1.5 2.0 

\\§ x - <]) II fn.n. x,y) 

FIG. 1. Probability distribution of \\<j> x — (f> y \\ for all n.n. 
x, y on a lattice of size L = 1024 at T = 0.5. In the inset: 
fraction of bonds that connect Tt + and Tt~ sites as a function 
of e. 

In Fig. H we show the distribution of sizes for the 
three types of clusters, S (circles) and H ± (triangles). 
We present the data for each cluster size C up to sizes 
C = 100; for larger values of C we show the results av- 
eraged over bins [InC — r;/2,lnC + 77/2] with r\ w 0.5. 
The quantity P(C) dC is proportional to the probability 
of finding a cluster with size between C and C + AC . 
The line with triangles is continuous and it ends at 
InC = 10.2. On the other hand the curve with circles, 
which represents the S clusters, displays two parts: a con- 
tinuous line ending at InC = 9.2 and a separated point 
at InC ps 13.1. This isolated point is the one determined 
by the percolating cluster. Although not explicitely de- 
picted, this point has an horizontal error bar which re- 
sults in the slight horizontal spreading of circles. This 
error bar is a remnant of the fact that on truly infinite 
lattices the size of the percolating cluster must be infinite. 
Notice that the continuous part of the curve with circles 
becomes steeper just before ending at InC = 9.2: this is 
another indication for the existence of a percolating clus- 
ter |pfj| because it means that all clusters beyond some 
size prefer to be absorbed by the percolating S cluster. 
This figure has been obtained by working on a lattice size 
L = 1024 and choosing ft — (1,0,0). Completely anal- 
ogous results are obtained with other sizes L and any 
other unit vector ft. We then conclude that there always 
exists a percolating cluster of type S for any versor ft. 
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In C 

FIG. 2. Size distribution of S (circles) and Tl ± (triangles) 
clusters on L — 1024 at T = 0.5 and e = 1 on a log-log scale. 
In the inset: ratio perimeter/size of clusters as a function of 
their size. 

In Fig. H an example of a type S percolating cluster 
taken from a single thermalized configuration is shown. 
The percolating sites are coloured in black. To render the 
fine details of the clusters as clear as possible we show 
only a piece of size 512 x 512 taken out from a configura- 
tion on a lattice with L — 1024. The figure corresponds 
to the percolating cluster found in the plane perpendic- 
ular to ft = (1, 0, 0). It is clear from this figure that the 
clusters present a high degree of roughness. 

The property of roughness is made explicit in the in- 
set of Fig. pi The ratio "perimeter of cluster" / "size of 
cluster", B/C, is displayed against InC for all 7i ± (tri- 
angles) and S (circles) clusters (irrespective of the fact 
that they percolate or not). For bidimensional compact 
objects B/C oc 1/yC when C — » 00. In our case this ra- 
tio tends to a constant (which is almost 1) when C — > 00 
and this is an indicative that our clusters present a rough 
border, suggesting that they have a fractal structure. 

In the inset of Fig. [2] we see that the ratio B/C goes 
asymptotically to a constant which seems to be a bit 
larger for the S clusters than for H . We think that this 
is because the S clusters surround the 7i ± clusters (see 
the introduction) and consequently they need to have an 
additional boundary. 

We have also studied the ratio among the size of the 
largest (always percolating) S cluster, Ms = maxjCg}, 
and the size of the largest 7i ± , Mu = max{C^±}, as a 
function of the lattice size L. Let us call R such ratio, 



If R increases with L for a fixed temperature, we can rea- 
sonably infer that the above described percolation prop- 
erties survive the thermodynamic limit. The data for R 
and the maximum size of S and Tt^ clusters are shown in 
Table 1 from the four lattice sizes we have used. They in- 
dicate that indeed the largest S cluster keeps percolating 
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when L increases. Notice that the percentage of lattice 
volume covered by the percolating cluster does not vary 
appreciably with the lattice size L. 




T > 0.5. 

We have also shown that, at least at T = 0.5, all clus- 
ters present a high degree of roughness recalling a fractal 
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It seems un- 



structure, (see inset in Fig. 
likely that, as suggested in |LC], such a dilute set of spins, 
almost lying on a plane, can enforce the system to behave 
like an effective 0(2) model (see pl|). 
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FIG. 3. Percolating S cluster (in black) from a thermalized 
configuration on a lattice of size L = 1024. 

TABLE 1. Ratio R and maximum sizes of clusters as 
a function of the lattice size L. 



L 


1024 


1250 


1550 


2048 


(M S ) /L* 


0.437(10) 


0.429(4) 


0.433(7) 


0.424(4) 


(M H ) /L* 


0.019(3) 


0.015(2) 


0.011(1) 


0.0072(7) 


R 


23(4) 


29(4) 


40(4) 


59(6) 



In conclusion, we have studied the percolation prop- 
erties of several clusters which can be defined on config- 
urations of the 2D Heisenberg model described by the 
hamiltonian in Eq. (|l|). Of particular interest are the 
clusters called S. To define them one has to introduce 
an arbitrary unit vector ft in the symmetry group space 
of the system, 0(3). When the spin cj> x at the site x is 
almost perpendicular to n, we say that this site belongs 
to some S cluster. The term "almost perpendicular" de- 
pends on the temperature T of the system (see above). 
For T = 0.5 and working on rather large lattice sizes 
L > 1024, we have given strong evidence that for every 
n one of these S clusters percolates on each thermalized 
configuration of the system (see Fig. ||). These percola- 
tion properties seem to survive the thermodynamic limit 
(see Table 1). 

This is an important conclusion because if such clus- 
ters do not percolate then one can prove Jll],[ll| that the 
system carries no massive particle, contradicting the ex- 
act calculation of the mass gap performed for this theory 
in Ref. M. Our results exclude this massless phase for 
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